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ABSTRACT 

We investigate the ability to discern between lognormal and powerlaw forms for the observed mass 
function of dense cores in star forming regions. After testing our fitting, goodness-of-fit, and model 
selection procedures on simulated data, we apply our analysis to 14 datasets from the literature. 
Whether the core mass function has a powerlaw tail or whether it follows a pure lognormal form 
cannot be distinguished from current data. From our simulations it is estimated that datasets from 
uniform surveys containing more than w 500 cores with a completeness limit below the peak of the 
mass distribution are needed to definitively discern between these two functional forms. We also 
conclude that the width of the core mass function may be more reliably estimated than the powerlaw 
index of the high mass tail and that the width may also be a more useful parameter in comparing 
with the stellar initial mass function to deduce the statistical evolution of dense cores into stars. 
Subject headings: stars: formation — ISM: clouds — ISM: structure 



1. BACKGROUND 

Stars form from dense cores of molecular gas and dust 
with sizes w 0.1 pc, masses between sa 0.5 and 50 



and d ensities > 10 cm 3 (see, e.g., Bergin & Tafalla 



20071) . Therefore the relationship between mass distri- 
butions of dense cores in star forming regions (core mass 
function, or CMF) and the stellar initial mass function 
(IMF) contains information regarding how observed sam- 
ples of cores evolve into stars. The observational simi- 
lar ity between the CM F and the IMF was first put forth 
by Motte et al. (19981, and since this time many other 
samples o t dense cores have been presented in this con- 

Enoch et al 



text (e.g 
et al. 2007 



Simp son et al.||2008| |Enoch et al |2008| |A lves 
Mutter fe Ward-Thompson|2007[ |Stanke et al. 



2006f|jon" nstone et al.||200l|). The qualitative similarity 



between the CMF and the IMF offers support for t he ac- 



cepted idea that stars form from dense cores (e.g. Lada 
et al. 2008| ). But what can we learn about the star tor 
mation process from this observed similarity? 

In a recent paper, we presented a series of numerical 
experi ments evolving distributi ons of cores into stellar 
IMFs ( |Swift fc Williams| |2008] . We find that a given 
CMF evolved according to different evolutionary path- 
ways produces variations in the resultant IMF that are 
insignificant in relation to the errors inherent in current 
samples of dense cores. Our results show that the form 
of the CMF in relation to the IMF indeed contains vital 
clues to the star formation process, but highlight the dif- 
ficulty in deducing how observed samples of cores evolve 
into stars based solely on the shapes of these distribu- 
tions. 

The central limit theorem applied to isothermal tur- 
bulence naturally predicts a lognormal probability distri 



bution in density (Larson 1973 Zinnecker 



|fc Fatuzzo 1996) that has been producec m computer 



simulations (e.g., [V azqucz-Semadeni 1994 



1984 Adams 



Nordlund fc 

|Padoan|19 99 ; Klcsscn 200 1|. Extensive surveys ot nearby 
star-forming regions have also shown a CMF shape con- 
sistent with a lognormal form (e.g., Enoch et al. 2008 
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Stanke et al][2006| . 

However, observed CMFs are typically characterized 
by one or more powerlaws. While the use of a powerlaw 
form to fit observed CMFs has never been rigorously jus- 
tified, this form is assumed based on its versatility and 
the e xpected similarity between the CMF and the IMF 



(e.g., Motte et al.[1998l. Recent data show that this ap- 



parent powerlaw behavior does not extend to very low 
masses but display s a turnover or break below a few Mq 
( Alves et al.||2007 l. Motivated by the observational sim- 
ilarity of the CMF and the IMF, theories have been re- 
cently developed describing physical scenarios that may 
produce a powerlaw distribution of molecular core masses 
that turns over towar d low-masses (IPadoan fc Nordlund| 
2 002||Shu et al.||2004| |Hennebelle fc Chabrier||2008[ ). 

To understand how, or if, dense molecular cores pro- 
duce the full spectrum of stellar masses, it is essential to 
understand the probability distribution function (PDF) 
from which the CMF is drawn. If the parent distribution 
has a powerlaw tail extending to high masses, this would 
support the idea that stellar mass is almost entirely de 



termined in the molecular cloud phase (e.g., Alves et al. 
2007 1 . If observed mass distributions of cores are found 



to be drawn from a purely lognormal PDF, the origin of 
the powerlaw distribution of stellar masses will remain an 
open question. A lognormal CMF would disfavor the idea 
th at massive stars form d irectly from massive cores (such 



as 



Krumholz et al. 2009), and may imply that massive 



stars form through mechanisms distinct from low-mass 



stars (e.g., Bonncll et al. 2004). Distinguishing between 



these two forms is complicated by the difficulty in mea- 
suring the CMF over large dynamic ranges and the fact 
that lognormal and powerlaw forms can look quite simi- 
lar over limited mass ranges. 

In 2006, Reid & Wilson concluded that observed CMFs 
are consistent with being drawn from the same par- 
ent distribution and that this parent distribution is 
consistent with the IMF. However, the results of their 
Kolmogorov-Smirnov (KS) tests applied to the median 
normalized CMFs suggest that the CMFs are indistin- 
guishable while their goodness-of-fit statistic based on 
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X 2 suggest that different core samples may have differ- 
ent preferred forms. Given these ambiguities and the fact 
that the KS test is somewhat compromised when applied 
to two distributions with normalized medians, we revisit 
the topic of discerning the form of the parent distribution 
to the CMF with a new approach. A renewed interest in 
this subject is motivated by new instruments coming on- 
line (e.g., Herschel, SCUBA2, ALMA) and surveys that 
will provide large, uniform samples of carefully selected 
cores 1 . 

In §[2] we outline our analysis strategy and then test 
our procedures on simulated data in §|3] The application 
of our procedures to 14 state-of-the-art datasets of dense 
cores in star-forming regions is presented in §[4j and we 
close with a discussion of these results in §[5]that includes 
suggestions for future observations. 

2. DATA ANALYSIS 

The goals of our data analysis are to (1) find the best 
fit parameters for lognormal and powerlaw models given 
the data, (2) assess whether or not a given model ade- 
quately describes the data, (3) select either the lognor- 
mal or powerlaw model as the preferred model, and (4) 
compute the confidence in that selection. The follow- 
ing sub-sections outline the details of each step used to 
achieve these goals. All procedures were written in the 
Interactive Data Language 2 . 

2.1. The Models 

The two models we choose to differentiate in our anal- 
ysis are motivated in §[l]and are described by probability 
distribution functions (PDFs) of powerlaw and lognormal 
form. The powerlaw PDF is given by 

Ppi(m) = C pl m _Q (1) 

for to > mo- The normalization constant is given by 

C p i = (a - l)™^ 1 

where too signifies the break or "turnover" below which 
the distribution does not follow powerlaw behavior. 
The lognormal PDF is given by 



Pin (m) 



Cln 
777 



exp 



(In 777 — fi) 2 . 
2a^ 



(2) 



with the normalization constant 

'lnmii m - /1 



Cln — 



;rfc 



aV2 



where the characteristic mass m c — and mn m is the 
minimum mass of the distribution, analogous to an ob- 
servational completeness limit. 

2.2. Fitting 

To avoid the dangers inherent in fitting data us- 
ing regression mo dels arising from, e.g., data binning 
( |Rosolowsky|2005| l, we use the method of maximum like- 
lihood to estimate the model parameters fo r both simu- 
lated and real data (see Clauset et al. 2007). 



1 e.g., The Gould's Belt Legacy Survey 
http: //www. iach .hawaii . edu/ JCMT/surveys/gb 



The best estimate of the powerlaw index for data with 
values above too is found analytically by maximizing the 
log- likelihood: 



Vln ( mi \ 



(3) 



where 77 is the number of data with 777 > m . Hatted 
quantities will signify the best fit parameters throughout 
the text . The value of too is also unknown and is found by 
minimizing the KS statistic betwee n the best fit model 
and the data as a function of too (Clauset et al. 112007 
§3-3). . 

The likelihood for the distribution in Equation can- 
not be maximized analytically. Therefore, we numeri- 
cally maximize the likelihood for the lognorm al model as 
a fun ction of a and \x using Powell's method (Press et al. 
20071. 



2.3. Goodness-of-fit 

Once the data are fit, we calculate the KS statistic 
between the best fit model and the data as a proxy for 
goodncss-of-fit. When one or more parameters of a model 
are determined from data, the probability value for the 
KS statistic, pks> can no longer be found directl y from 
the KS probability distribution (Lilliefors 1967). We 
therefore calculate pks by relating the KS statistic to 
the distribution of valu es generated from 5000 parametric 



bootstrap realizations (P ress et al. | 2007 Babu &; Feigel- 



son 



sstrap realizations ( rress 1 
2006] |Clauset et al ||2007| 



he value ot pxs represents the probability of getting 
a fit worse than the original given the model and the 
PDF of the distribution. Low values of pks, typically 
below 0.1 or 0.05, are a flag for a given model being an 
inadequate description of the data. 

2.4. Model Selection 

There are several approaches to selecting between mod- 
els. Bayesian model selection compares the evidence of 
each model computed as an integral over all free param- 
eters of the probability density, i.e., the likelihood func- 
tion. One of our free parameters, m , determines the 
number of data points to be considered by the models 
thus adding a difficult complexity to this approach. In- 
stead we use the likelihood ratio test, a similar approach 
that compares the peak probability densities of the mod- 
els. Taking the log of this ratio we obtain the parameter 



TZ = ^2 [In p p i (rrii) - \TLp\ n (m % 



(4) 



1DL; http : //www. ittvis . com/ProductServices/IDL . aspx 



where positive values of 1Z favor the powerlaw model and 
negative values favor the lognormal model. 

Our two models each have two parameters to and 
a for the powerlaw model, and /j, and a for the lognor- 
mal model. Therefore we do not expect there to be bi- 
ases introduced into this method due to different levels 
of complexity in the models. An additional benefit to 
using this test is that the significance of sgn(7?.) can be 
assessed by generating a distribution of 1Z values through 
bootstrap resampling of the data (for real data) or the 
PDF (for simulated data) . The fraction of 1Z values hav- 
ing the same sign as the original value is designated C-jz, 
the confidence level of the likelihood ratio. 
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3. APPLICATION TO SIMULATED DATA 

To ascertain how well our procedures discriminate be- 
tween powerlaw and lognormal distributions, we have 
carried out our analysis on a large number of simulated 
datasets. These simulations quantify how sensitively the 
model selection process depends on the sample size and 
mu m . 

Distributions with Nc samples above a completeness 
limit of m lim for A^ c = 50 x 2™ (n = 0,1, 2,. ..,8) and 
rnii m = 0.1, 0.5, and 1.0 Mq are drawn from the PDFs 
of § |2.1| For powerlaw data, we draw from the distribu- 
tion in Equation [l] for masses above mo while for masses 
between mn m and mo we draw from Equation |2] The 
powerlaw data are drawn such that the parent distribu- 
tion is smooth and continuous at uiq = /i — er 2 (l — a). 
Lognormal data are drawn from the distribution in Equa- 
tion^ For each sample size, N c , and completeness limit, 
WUnu 5000 distributions are realized and fit according to 
the procedures outlined in §[2j 

We use values a = 1.0 and fi = for all simulated 
data, consistent with CMFs measured in nearby star - 
forming regions ( Enoch et al.||2008 Stanke et al.| [2006). 
For powerlaw data we use a = 2.5, close to the Salpeter 
value of 2.35 (Salpeter 1955ft and also consist ent with 
data from the literature ( Keid fc Wilson]|2006b also see 
Table 

3.1. Likelihood Ratio Test 

Figure [T] displays C-r as a function of Nq given differ- 
ent values of mn m . The solid lines in the figure represent 
the confidence with which the underlying powerlaw form 
to the PDF can be discerned from a lognormal form. 
The confidence is an increasing function of Nq and a de- 
creasing function of mii m expressed here as the fraction 
of the characteristic mass of the distribution, my lm /m c . 
The dependency on Nc is expected and reflects the in- 
creasing amount of information available to the likeli- 
hood ratio test. The dependency on mi; m arises because 
a lognormal form can closely approach a powerlaw form 
over a limited range for large negative values of fi and 
large positive values of a. Therefore, as more data are 
included below uiq stronger constraints can be placed on 
the values of /t and a thereby creating a larger discrep- 
ancy between a lognormal model and the high mass tail 
of the simulated data. We find that a lognormal form to 
the CMF can be ruled out at the 95% confidence level 
for a survey of intermediate sensitivity and sample sizes 
greater than w 500. 

The dashed lines in Figure [l] represent the confidence 
level with which a pure lognormal form to the PDF can 
be distinguished from a one with a powerlaw tail. We 
see that these curves are monotonically increasing for 
Nc > 100. The upturn in the lognormal data curves 
at low Nc is due to the effects of fitting a powerlaw to 
a small number of data points. As Nc increases, so do 
the best fit values, mo and a, thereby making a distinc- 
tion between models more difficult than for the case of 
powerlaw data. If mo is forced to lie at some fraction 
of m c , then results similar to the a powerlaw data (solid 
lines) are obtained. The dependency of C-r on mi; m is 
weaker for the lognormal data than the powerlaw data 
since lower values of mi; m produce only slightly better 
lognormal fits. 




Powerlaw Data 
Lognormal Data 



m lim /m c = 0.1 
m lim /m c = 0.5 
m,| m /m c = 1.0 



100 1000 10000 

Fig. 1. — The confidence with which the likelihood ratio test can 
discern between a powerlaw and lognormal model for data gener- 
ated from a powerlaw (solid lines) and lognormal (dashed lines) 
probability distribution function as a function of the sample size, 
Nq. Data with completeness limits of 0.1, 0.5, and 1.0 times the 
characteristic mass are shown in black, red, and blue, respectively. 



3.2. Goodness- of- fit Test 

Figure [2]shows the average value of pks as a function of 
Nc generated by fitting the incorrect model to the simu- 
lated data. For both the powerlaw and lognormal data, 
the incorrect model produces poorer fits as the sample 
size increases. Similar to the likelihood ratio test, it is 
more difficult to rule out a powerlaw model from lognor- 
mal data. However, the discrepancy between the critical 
values of Nc is smaller using this test. 

For powerlaw data, (pks) is a decreasing function of 
mii m . This dependency comes about since the KS test 
is most sensitive to the median of a distribution and the 
powerlaw and lognormal models have median values that 
become closer as mi im decreases. For lognormal data, 
(j>ks) is not a monotonic function of my lm for all values 
of N c . 



Powerlaw Data 
Lognormal Data 




Fig. 2. — The average KS probability resulting from fitting a log- 
normal model to powerlaw data (solid lines) or a powerlaw model 
to lognormal data (dashed lines). Color conventions are the same 
as in Figure [l] 
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It may be expected that the curves of (pks) should 
approach 0.5 as Nq — > 0, meaning that the two models 
are indistinguishable for small sample sizes. However, 
it can be seen that the curves representing lognormal 
data (dashed lines) overshoot this value for Nq between 
about 50 and 200. As in § |3.1| this results from fitting 
a powcrlaw to a small number of samples and represents 
a bias for lognormal data to look marginally more like a 
powerlaw for Nq in this range according to our goodness- 
of-fit test. 

From these simulations it can be concluded that fit- 
ting a powerlaw model to lognormal data or conversely a 
lognormal model to powerlaw data will consistently yield 
very low values of j?ks for Nc greater than several hun- 
dred to a thousand. 

4. APPLICATION TO LITERATURE DATA 

Table [l] lists the 14 datasets to which we apply the 
analysis of §|2] along with their descriptive parameters. 
Various techniques are used to observationally identify 
the cores in these datasets. The most common techniqu e 
is the use of threshold contou ring (e.g., Alves et al.|2007|), 
thoug h gaussian fitting (e.g., Nutter & Ward-Tho mpson 
2007 1 and wavelet decomposition (e.g., Motte et al.|1998| ) 
are also used. 

The results of our fitting procedures, goodness-of-fit 
tests, and model selection techniques are displayed in 
Table [2j For a majority of the datasets, neither a lognor- 
mal form nor a powerlaw provide a significantly better 
description of the data. Five datasets out of 14 have val- 
ues of pks and C-jz that are potentially significant. These 
values are underlined in the table for emphasis. 

Datasets 3, 11 and 14 disfavor a powerlaw model due 
to a very low pxs value while showing various levels of 
support for a lognormal model in the likelihood ratio. 
Dataset 12 disfavors a lognormal form in the pks statistic 
while showing moderate support for a powerlaw form in 
the likelihood ratio test. Dataset 6 is not adequately fit 
by either model. 



TABLE 1 
Literature Data 



5. DISCUSSION 

5.1. Does the CMF Have a Powerlaw Tail? 

Ambiguous results for more than 60% of the datasets 
and conflicting results for the remaining ~ 40% empha- 
size the difficulty in discerning the form of the CMF. 
Given the results in Table [2] we conclude that a pure log- 
normal PDF and a PDF with a powerlaw tail cannot be 
distinguished from current core datasets. Furthermore, 
neither the existence of a powerlaw tail to the CMF nor 
the uniformity of the CMF from region to region are ev- 
ident. 

It is not clear to what degree systematic differences in 
each dataset — differing observational techniques, phys- 
ical resolution, filtering, core identification algorithms, 
and parameters used in deriving core mass (e.g., k v ) — 
contribute to the variations seen from dataset to dataset, 
and it is possible that a uniform re-analysis of the 
datasets could reveal a more definitive relationship be- 
tween the shapes of the CMFs. Also, the free parameters 



Dataset ID 


Region 


Wtot 


miim 


N c 


Ref. 


1 


Pipe JNcbula 


loa 


i n 


— Q 


i 

1 




Pipe Nebula 


1.1 1 


U.D 




o 


o 


Bolocam Cores (SL) 


lUo 


U.o 




■3 

O 


4 


Bolocam Cores (PS) 


92 


0.8 


15 


3 


5 


Bolocam Cores (All) 


200 


0.8 


114 


3 


6 


p Ophiuchus 


55 


0.1 


40 


1 


7 


p Ophiuchus 


143 


0.1 


120 


5 


8 


p Ophiuchus 


62 


0.1 


51 


6 


9 


Orion A 


172 


0.3 


170 


7 


10 


Orion 


203 


1.0 


149 


8,9,10 


11 


M17 


121 


3.0 


108 


11 


12 


NGC 7538 


77 


15 


15 


12 


13 


Cygnus X 


129 


5.3 


125 


13 


14 


NGC 6334 


181 


32 


131 


11 



(1) |Alves et al.| ||2007|l. (2) |Rathborne et al.| 
al.| ^UUalrT^ ) Johnsto ne et al.| pJUo), (5) 
"et al.N2UU6 |TT ynCToEte et al.TfLyiJ81, ('V') [Nutferir'Ward- | 

Thompson (2007), (8) Johnstone et al. (2001), (9;johnstonc et al. 
(20U6), (10) Johnstone & Bally (2006). (11) Hcid & Wilson (2006a), 
(12) Rcid & Wilson (2005k, (13) llVlotte et al.N2007|, (14) |Munoz| 

|itlr.||2uuyt 



References. 
( 200 9), (3) lEnoch e t 
Stankc 



§[3] The literature datasets of Table [l] have 
vrhilc the results of S[3] suggest that more than 



are unconstrained in our fitting procedures, and in some 
cases this leads to fits that could be seen as unreasonable 
(e.g., the low values of jj for datasets 8 and 9). Judicious 
constraints on the fitting parameters could therefore lead 
to more decisive results. 

However, our conclusions are consistent with the re- 
sults from 

N c < 150 wtme the results of §[3]s 

500 samples are needed to dermitively discern between 
these two models. These results highlight the need for 
larger, uniform survey data to produce a statistically sig- 
nificant result regarding the functional form of the CMF. 

5.2. Comparing the CMF to the IMF 

The utility in analyzing the form of the CMF is in 
its relation to the stellar IMF to address the question of 
how the final mass of a star is determined. The similarity 
of the CMF to the IMF has inspired the idea that the 
masses of stars are determined via the fragmentation of 
dense molecular gas in the Galaxy. It seems as though 
this must be true at least in part since it is known that 
stars form in dense, self-shielded molecular gas. However, 
whether or not the mass distributions of molecular cores 
observed in star-forming regions will accurately represent 
the mass distribution of stars to form from them is still 
debatable. 



It was shown by Swift & Williams ( 2008[ ) that different 
evolutionary pathways from cores to stars produce vari- 
ations in the form of the resultant IMF. The powerlaw 
slope was shown to be quite robust except under the most 
extreme evolutionary scenarios due to the scale-free na- 
ture of powerlaw distributions. The width however was 
shown to be a more sensitive indicator of core evolution. 
Given the small dynamic range over which the CMF can 
be measured in nearby star-forming regions, we empha- 
size the importance of obtaining a sound measurement 
of g in survey data as this parameter may be more pow- 
erful than a in constraining how observed distributions 
of cores evolve into stars. 
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TABLE 2 

Literature Analysis Results 



Powcrlaw Fit 



Datasct ID 



mo 



1 2 9 +0 ' 6 

2 3.ot°;g 

n o 7+0.3 

z -'-0.3 

4 3 4+ ' 5 

5 4.HJ;° 

6 i.t;;;:' 

7 2.8t°;| 

8 2.0±q 2 

9 i.9±8ii 

10 2.7±S"| 

12 2 1+ - 4 

13 2.lto 2 

14 i-rtft 



2.5 
2.4 
0.9 
2.5 
2.7 
0.2 
0.6 
0.1 
6.4 
3.1 



+0.04 

1.1 
+0.3 
1.3 
+0.7 
0.1 
+0.1 
0.9 
+0.1 
1.5 
+0.5 

0.01 
+0.04 

0.2 
+0.1 
-0.02 
+ 1.8 
5.3 
2.4 

o.e 



14 



+6 



68 
21 
47 



+4 

45 



PR'S Mail 3 



Lognormal Fit 



5 

+41 
13 



0.75 
0.43 
0.00 
0.72 
0.32 
0.00 
0.85 
0.15 
0.67 
0.21 
0.00 
0.57 
0.66 
0.09 



27 
30 
65 
22 
28 
37 
43 
46 
37 
69 
59 
22 
65 
111 



PA'S 



+0.6 
-2.6 

2+ ' 3 
u -^-0.6 

0+ ' 3 

7+ ' 2 
,+0.2 
-0.5 
-,+0.2 
-0.3 
1+0.2 
-0.3 
,+2.1 
-2.0 
4+1.4 
-3.9 
7+0.2 
-0.4 
q+0.2 
-0.3 

2.7+ 1 ' 1 



-0.1 



0.2: 

-0.9 1 
-1.1" 1 
-4.lj 
-5.4j 
0.7j 
2.6^ 



2.7 

2.6 



3.0 

+0.3 
0.5 
+ 1.3 
3.9 



1+0.8 
-0.3 
h0.3 
-0.2 
7+0.3 
-0.2 
7+0.3 
-0.1 
3+0.2 
-0.1 
(-0.3 
-0.2 
{-0.2 
-0.1 
h0.8 
-0.7 
hl.l 
-0.5 
1+0.2 
-0.1 

1 1+ - 2 
J — L -0.2 
hi. 3 
-0.6 
j+0.3 
-0.2 
hi. 2 
-0.5 



l.oj 

0.7j 
0.7j 
0.8j 
0.9^ 
1.0 4 
1.9j 
3.4j 
1.0- 1 



1.9 -1 
1.3j 
2.2 4 



miim/m-c H Cn 



n on 


i . i 


1 L 


tv/o 


U.ol 


U.u 


1 L 


(1/0 


0.42 


0.8 


LN 


91 70 


0.68 


0.4 


PL 


68% 


0.54 


0.7 


PL 


60% 


0.06 


0.2 


LN 


77% 


0.45 


0.3 


PL 


85% 


0.20 


5.9 


LN 


69% 


0.94 


200 


PL 


56% 


0.45 


0.5 


PL 


68% 


0.59 


0.2 


LN 


90% 


0.07 


1.0 


PL 


84% 


0.14 


2.2 


PL 


83% 


0.74 


2.3 


LN 


74% 



a Number of cores with m > mo. 



The width of the IMF has been measured to be be- 



tween 0.3 and 0.7 dex (Mil ler fc Scalo ||1979| |Chabrier 
2003 Bochanski et aL 2009). The measured values of a 



in Table [2] tend to be mgher than this, approximately be- 
tween 0.7 and 2. If this were confirmed to be the case in a 
more sensitive, uniform survey, this would indicate that 
an additional mass preference occurs in later stages of 
gravitational collapse likely through fragmentation into 
several stars at or near the local Jeans mass. 

5.3. Bettering Our Understanding of the CMF 

It is a difficult task to automate the identification 
of large numbers of pre-stellar molecular cores in star- 
forming regions, and there is no consensus on how to 
define a core observationally. H owever, methods are cur 



rently bein g improved (Kainulainen et al. 2009 Pineda 



et al.|2009ft and new stra tegies are being developed (e.g., 
olo- 



Rosolowsky et al.||2008[ ). It has been rece ntly shown 
that the phy sical "characteristics of cores (e.ff., |Rathborne| 
et al. | 2009|), the p otential variability of core lifetimes 
( jciark eF al. 2007) as well as the effects of statistical 
errors (e.ff., |Koen &; Kon dlo 2009 Rosolowsky|2005j ) are 
also important to consider when interpreting observed 
CMFs. 

From this work, we present three new suggestions 
aimed to help maximize the utility of future surveys of 



dense cores in star-forming regions. One, is that sam- 
ple sizes of greater than 500 above the completeness 
limit should be sought so that a lognormal and powerlaw 
tail to the CMF can be discerned. The second sugges- 
tion is that the survey should be as uniform as possible 
such that systematic errors can be minimized. The third 
suggestion is that a sensitivity limit of better than one 
half of the peak mass of the distribution is desirable such 
that the width of the CMF can be reliably determined 
and compared to the width of the IMF. 

A survey such as the Gould's Belt Legacy Survey to be 
conducted with SCUBA2 on the James Clerk Maxwell 
Telescope may very well reach these goals, and the value 
of a is likely to be measured with unprecedented accu- 
racy. However, this survey will only cover low to in- 
termediate mass star-forming regions and it may not be 
possible to determine a robust value of a. For this, we 
may need to wait for the exquisite resolution and sensi- 
tivity of ALMA to probe nascent regions of massive star 
formation. 
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